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Abstract 

Motivated by the recent computations of the quasinormal 
frequencies of higher dimensional black holes we exactly 
calculate the quasinormal frequencies of the Dirac field 
propagating in D-dimensional (D > 4) massless topologi- 
cal black hole. From the exact values of the quasinormal 
frequencies for the fermion and boson fields we discuss 
whether the recently proposed bound on the relaxation 
time of a perturbed thermodynamical system is satisfied 
in D-dimensional massless topological black hole. Also we 
study the consequences of these results. 

Keywords: Quasinormal modes, massless topological 
black hole, Dirac field, relaxation time 

Motivados por el calculo de las frecuencias cuasinor- 
males de agujeros negros cuyo mimero de dimensiones 
D es mayor o igual a cuatro, en el presente arti'culo 
calculamos exactamente las frecuencias cuasinormales del 
campo de Dirac moviendose en el agujero negro topologico 
de masa cero con D > 4. Usando los valores exactos de 
las frecuencias cuasinormales para los fermiones y bosones 
discutimos si el limite, recientemente propuesto, sobre 
el tiempo de relajamiento de un sistema termodinamico 
perturbado se satisface en el agujero negro topologico de 
masa cero con D > 4. Adicionalmente estudiamos algunas 
consecuencias de estos resultados. 

Descriptores: Modos cuasinormales, agujero negro 
topologico de masa cero, campo de Dirac, tiempo de rela- 
jamiento 

PACS: 04.70.Bw, 04.50. Gh, 04.70.Dy 

1 Introduction 

The physical systems for which we exactly solve their equa- 
tions of motion can be expected to play a significant role 
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in several research lines. For these physical systems we 
exactly calculate the physical quantities that for other sys- 
tems we calculate by using approximate methods. Also in 
many research areas the physical insight that is obtained 
by studying the exactly solvable systems can be used to in- 
fer some details about the behavior of more complex phys- 
ical systems. 

The quasinormal modes (QNM) of a black hole are so- 
lutions to the equations of motion for a classical field that 
satisfy the appropriate radiation boundary conditions at 
the horizon and at the asymptotic region. The quasinor- 
mal frequencies (QNF) of a field are valuable quantities 
since these are determined by a few parameters of the black 
hole and the field Q]-[3], for example, the QNF of the Kerr- 
Newman black hole are determined by the mass, angular 
momentum, and charge of the black hole and the mode of 
the field. Hence if we measure the QNF of a field then we 
can infer the values of the mass, angular momentum, and 
charge of the Kerr-Newman black hole. 

Also the QNM allow us to study the linear stability of 
the black holes, because if we find QNM whose amplitude 
increases in time, then the black hole may be unstable [I]- 
[3]. Recently the QNM have found applications in several 
research lines. For example, a) the AdS/CFT correspon- 
dence of string theory [2], [I], 0, b) the determination of 
the area quantum of the black hole event horizon [6] , [7] , c) 
the expansion of functional determinants in some thermal 
spacetimes [8] , [9] , d) the expansion of the "distant past" 
Green functions used in self- force calculations [IT)] . 

For many relevant spacetimes their QNF must be calcu- 
lated approximately, hence we use numerical methods or 
perturbation methods Q~]-[3]- Nevertheless, recently ex- 
act calculations of the QNF for several spacetimes have 
been presented. Among these we enumerate the follow- 
ing, a) three-dimensional static and rotating BTZ black 
holes [5], [H]-[I3], b) three-dimensional charged and ro- 
tating black holes of the Einstein-Maxwell-dilaton with 
cosmological constant theory [2]-[T7], c) two-dimensional 
dilatonic black hole [18] , [19] , d) five-dimensional dilatonic 
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black hole [18] , [19] , e) D-dimensional de Sitter spacetime 
(D > 3) EU]-[26], f) BTZ black string [27], g) Nariai space- 
time [28 o In the following paragraphs we comment on 
another Z?-dimensional anti-de Sitter black hole for which 
the exact values of its QNF have been calculated. 

We notice that the AdS/CFT correspondence of string 
theory motivated many studies on the QNF of anti-de Sit- 
ter black holes 0, [I], [5], because this correspondence 
proposes that the QNF of the anti-de Sitter black holes 
determine the relaxation time of the dual conformal field 
theory [1] , [5] . See Ref. [5] for an explicit verification of this 
proposal in three-dimensional rotating BTZ black hole. 

Furthermore we recall that in asymptotically anti-de Sit- 
ter spacetimes there are solutions of the Einstein equations 
that represent black holes whose horizons are negative cur- 
vature Einstein manifolds [30]-[36]. These solutions are 
usually known as topological black holes and for some of 
these solutions the mass parameter can assume negative 
or zero values [30]-[36], 

Among these exact solutions of the Einstein equations 
there is one that has attracted a lot of attention. It is 
the asymptotically anti-de Sitter black hole whose mass is 
equal to zero (30]-[3J)]. In the rest of the present paper 
we call it the massless topological black hole (MTBH). 
According to Ref. [37], we can consider the MTBH as a 
higher dimensional generalization of the three-dimensional 
static BTZ black hole and we expect that it will play a 
significant role in future research. 

The metric of the £>-dimensional MTBH is simple and 
as a consequence many of its physical properties can be 
calculated exactly [3_7]-[4T]. For example, the QNF of the 
gravitational, Klein Gordon, and electromagnetic pertur- 
bations were calculated exactly in Refs. [37], [38], and Sec. 
6 of Ref. [39], respectively. Also its stability against the 
three types of gravitational perturbations was proven in 
Refs. [40], [41]. For numerical and analytical computa- 
tions of the QNF for other topological black holes see Refs. 
@2]-[48]. 

Here we exactly calculate the QNF of the Dirac field 
evolving in D-dimensional MTBH and thus we extend the 
results of Refs. [37]-[39]. The computation of the QNF 
for this fermion field is interesting because in some back- 
grounds the Dirac field behaves in a different way that the 
boson fields, for example, it is well known that in a rotat- 
ing black hole the Dirac field does not show superradiant 
scattering [49]-[52] in contrast to boson fields [53] ■ Also 
notice that the QNF of the Dirac field allow us to discuss 
some additional details about the behavior of the MTBH 
under perturbations. 

Note that in higher dimensional spacetimes, for the 



x Wc notice that in Ref. |29| Saavedra presented an exact expres- 
sion for the QNF of Unruh's acoustic black hole. The expression 
used in that reference for the effective metric of Unruh's acoustic 
black hole is valid near the horizon. For the asymptotic region of 
Unruh's acoustic black hole it is probable that we need to use a dif- 
ferent approximation of the effective metric. Thus we believe that 
this problem deserves additional study. This issue was pointed out 
to the Author by the Referee. 



Dirac field we only know the QNF reported in Refs. [19] , 
[26] , [54] -[57], thus for this fermion field its resonances 
have not been studied as extensively as for other fields. 
Hence this paper extends our knowledge on the QNM of 
the Dirac field in higher dimensional black holes. 

This paper is organized as follows. In Sec. [2] we find 
exact solutions to the Dirac equation in D-dimensional 
MTBH and using these solutions we exactly calculate the 
QNF of the Dirac field. Exploiting these results we enu- 
merate some facts about the behavior of the MTBH under 
perturbations. In Sec. [3] we investigate whether the fun- 
damental QNF of the MTBH satisfy the bound recently 
proposed by Hod in Ref. [58]. In Sec. H] following Chan- 
drasekhar [53], in MTBH we write the Dirac equation as 
a pair of Schrodinger type differential equations and iden- 
tify the effective potentials. Finally in Sec. [5] we discuss 
the results obtained. 



2 QNF of the Dirac field 

The line element of a Go-2-symmetric spacetime may be 
written as |59 



ds 2 



F(r) 2 dt 2 - G{r) 2 dr 2 - H(r) 2 dE 2 D _ 2 



(1) 



where F(r), G(r), and H(r) are functions only of the coor- 
dinate r and dY? D2 denotes the line element of a (D — 2)- 
dimensional Gd- i -invariant base spacetime Y>d-2, which 
depends only on the coordinates <f>i, i — 1,2,. . . ,D — 2. 

Our aim is to calculate exactly the QNF of the Dirac 
field evolving in Z?-dimensional MTBH. Thus first we ex- 
plicitly write the Dirac equation 



iflip = mip 



(2) 



in MTBH to find its exact solutions. Note that we follow 
the usual conventions, thus in formula ([2]) the symbol JC 7 
denotes the Dirac operator, m stands for the mass of the 
Dirac field, and tp denotes the spinor of dimension 2^ D / 2 \ 
where [D/2] denotes the integer part of D/2 [56], [60]-[65], 
As is well known, in a D-dimensional Gi)-_2-symmetric 
spacetime with line element (fTJ , the Dirac equation reduces 
to a pair of coupled partial differential equations in two 
variables (see for example Eqs. (30) of [53] and Refs. [6"0]- 
[65] for more details) 

F ( F \ 

df4>2 ~ Q d r^2 = ~ lmF J V'l, 

F ( F \ 

d t ipi + Q d r^i = - I + imF\ -02, (3) 

where k stands for the eigenvalues of the Dirac operator on 
the manifold Sd-2 with line element dY? D _ 2 and the func- 
tions "01 and ?/>2 are the components of a two-dimensional 
spinor ip2D which depends only on the coordinates (t, r) of 
the Gc_2 _ symmetric spacetime with line element (fTJ), that 



i>2D(r, t) 



Mr,t) 



(4) 
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We point out that in Eqs. ([3]) and in the rest of this paper 
we write the functions ipi(r,t), ip 2 (r, t), F(r), G(r), and 
H(r) simply as ipi,ip2, F, G, and H, respectively. We shall 
use a similar convention for the functions to be denned in 
the rest of the present work. 

The line element of the Z?-dimensional MTBH is given 
by [3Q]-[36] 

r 2 \ , „ dr 2 



ds z 



1 



r 2 dX 2 D _ 2 , (5) 



where r € (L, 
stant A by 



-oo), L is related to the cosmological con 
(D-1)(D -2) 



L z = -- 



2A 



(6) 
2> 



and dY? D 2 stands for the line element of a (D 
dimensional compact space of negative curvature 
[2Q]-[36]. Notice that the (t,r) sector of the line ele- 
ment ([5]) for the MTBH is similar to that of the three- 
dimensional static BTZ black hole with mass M = 1. 
Taking into account this fact, it was proposed that the 
Z?-dimensional MTBH §5§ is a higher dimensional gener- 
alization of the three-dimensional static BTZ black hole 
[37]. 

The QNM of the MTBH are solutions to the equations 
of motion for a field that are purely ingoing near the event 
horizon and since this black hole is asymptotically anti-de 
Sitter, we impose that at infinity the radial functions go 
to zero (Dirichlet's boundary condition) [37], [39]. In this 
section we compute the QNF of the Dirac field propagating 
in D-dimensional MTBH to find out about the behavior of 
this black hole under fermion perturbations and compare 
with its behavior under boson perturbations. We note that 
the results of this section are an extension of those already 
published in Refs. [37]-[39]. 

The line element of the D-dimensional MTBH ([SJ 
has the same form that the line element of the Go -2- 
symmetric spacetime ([l}. Thus making the appropriate 
identifications we get that the functions F, G, and H for 
the MTBH are equal to 

1 / r 2 \ 1/2 
F = - = (-l + T2 ) , H = , (7) 

Therefore in D-dimensional MTBH the coupled partial dif- 
ferential equations ([3]) reduce to 

3*V>2 - £j^d s ih = (z 2 I) 1 ' 2 (i£ - m) ifc, 

fl^i + '—^-d^! = -(z 2 - 1) 1/2 (j^ + im) fh, (8) 

where z = r/L and therefore z 6 (1, +oo). In what follows 
we write in detail the procedure used to solve exactly Eqs. 
©• 

Choosing for the components ipi and i\) 2 a harmonic time 
dependence of the form 

1 p 1 (z,t)=R 1 (z)e- iut , 

M*>*) = R2(z)e-^, (9) 



and defining u) = u>L, rh = mL, and K = —in we get that 
the system of partial differential equations ((H) transforms 
into the coupled system of ordinary differential equations 
for the functions i?2 and R\ = —iR\ 

(z 2 - l)^ 2 - + iuR 2 = (z 2 - l) 1 / 2 (*J - m) Rt, (10) 

(z 2 1)^1 - iQR 1 = -(z 2 I) 1 ' 2 + fhj R 2 . 

If we make the following ansatz for the functions i?i 
and R 2 (see formulas (26) of Ref. [26] for a similar ansatz 
for the radial functions of the Dirac field evolving in D- 
dimcnsional de Sitter spacetime) 

R l {z) = {z 2 -l)- 1 l\z + l) 1 l 2 R l {z), 

R 2 (z) = (z 2 - l)-^(z - l) 1 ' 2 R 2 {z), (11) 

then we find that the functions R\ and R 2 satisfy 

, o , di?2 i i •, ~ / iK \ , , ~ 
(z 2 - l)-^- + (iw + \) R 2 = ( — - m J (z + l)R u 

(z 2 -l)^-(iQ+i)R 1 = -(ij-+ fhj (z - 1)R 2 . 

(12) 

Next, we define the functions fi and f 2 by 

h(z) = Rt(z) + R 2 (z), f 2 {z)=R l {z)-R 2 {z), (13) 

to obtain that these functions must be solutions to the 
coupled system of ordinary differential equations 

(z 2 1)^ + [mz h = (iu+±+iK-m) f 2 , 

( z2 - l )^-{^ z - l -^)h = (^ + \- l K + m)f 1 . 

(14) 

From Eqs. (TT4"]) we obtain that the functions /i and f 2 
satisfy the decoupled ordinary differential equations 



(z 2 - if^fy + 2z(z 2 - + (** - 1) I m + 



df 



1 , (Ji 



fh z — 2miK / : 



iK 



(^+i) 2 -^-m) 2 )/ 1 , 



dz 



iK 
1 



z 



m 2 z 2 - 2miK - ) /: 



K 



(^+i) 2 -{iK-rh) 2 )f 2 



(15) 



To solve Eqs. ([T5]) we make the changes of variables x = 
z 2 and u = (x — l)/x, and take the functions f\ and f 2 in 
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the form 



where 



f 1 (u) = u Bl (l-u) Fl R 1 (u), 
f 2 (u) = u B *(l-u) F *R 2 (u), 



B\ —B 2 



Fx 



Fo 



iui I 1 
2^4- 



iui 1_ 

2 A' 



Now we use these results to compute the QNF of the 
Dirac field exactly. First let us study the function fx. We 
choose the quantities Ci, B\, and F\ as C\ = 1/2 + iK/2, 
(16) B x = iti/2 + 1/4, and Fi = 1/4+ \/m 2 - m + 1/4/2. If 
we assume that the quantity c\ is not an integer, then we 
obtain that the function fx is equal to 



(17) 



3A /m 2 -m+ 3 , 



L lm 2 +m+ j, 



to find that the functions R\ and R 2 must be solutions of 
the hypergeometric differential equation [66) . [6 7) 

u(l - u)^4 + (c - (a + b + l)u)^- - abf = 0. (18) 



du 



If the parameter c is not an integer, then the solutions of 
Eq. (|18[) are given in terms of the standard hypergeometric 
functions 2 Fi(a, b; c; it) |BB], [57] . 

For the functions i?i and i?2 the quantities a, b, and c 
of Eq. (|18|) are equal to (a,, &i, and c, correspond to the 
function Ri, i = 1,2) 



«i 

h 
ci 

a2 

C2 



B\ + Ci 

Si - Cr 
2Bi + 1, 
B 2 + C 2 

B 2 — C 2 
2B 2 + 1, 



rh 2 — rh + 



rh 2 — m + 



(19) 



rh 2 + m + j , 



where the quantities C\ and C 2 take the values 

2 ~ 2 ' 
Ci = ^ C 2 = 



2 ' 



2 ' 



iK 
2 ' 



(20) 



At this point we notice that the coordinate x lies in the 
range x € (l,+oo). Hence the variable u satisfies u € 
(0,1). Also the tortoise coordinate of the MTBH is [33] 

r, = / I — 1 + — | ^dr = — L arccoth(2), ( 21 ) 



r 

I 2 



thus r* e (— oo,0), r* — > — 00 near the event horizon and 
r„ — > near infinity. From these definitions of the coordi- 
nates u and r* we get that 

3 2r„/L 



as 
as 



—00, 
0, 



u « e 

ttfSl. 



and 



(22) 



/1 = (1 -«) Fl {Di^ /2+1/4 2 ^i(ai,&i;ci; M ) (23) 
+ Eiu- ia, / 2 - 1 / 4 2 Fi (01 - ci + 1, bi - a + 1; 2 - c i; tt)} , 

where Bi and Ei are constants. Taking into account ex- 
pressions (|22p we find that near the horizon the function 
/1 behaves as 



fi 



he 



iuir„+r,/( 2 L) 



+ Eie 



-itor* — r* / (2L) , 



(24) 



thus in order to have a purely ingoing wave near the event 
horizon we must impose the condition Bi = 37 -[39 . 
Hence the function fx becomes 



x ( fl i - ci + lj&i-ci + l;2-ci;u) (25) 

= EiM" l ' : ' /2 " 1/4 (l - u ) 1 /4+V™ 2 -™+ 1 /4/2 

x 2-F 1 i(ai,/3i;7i;?/). 

We recall that if the quantity c — a — b is not an inte- 
ger then the hypergeometric function 2 Fx(a,b;c]u) satis- 
fies [661, EH 



2 Fi(a, 6;c; u) 



_ r(c)r(c - a - b) 
~ r(c - a)r(c - b) 
x 2-F1 (a, b; a + b + 1 — c; 1 — u) 
T(c)T{a + b - c) 



-(1-u) 



c—a — b 



(26) 



r(a)r(6) 

x 2-Fi (c — a, c — fe; c + 1 — a — 6; 1 — u) , 

where T(a;) stands for the gamma function. Hence if the 
quantity 71 — a% — fix is not an integer, then we write the 
function fx of formula (|25|) as 



fx = Exu 



r( 7l )r( 7l - ax - px) 



r( 7 i - ai)r( 7 i - Px) 



(27) 



X (1 - u )l/4+V™ 2 -» l + 1 /4/2 

x 2 Fi(ai,/3i;ai + ^1 + 1 - 71; 1 - «) 

, r( 7l )r(ai + /3 1 - 7l ) 



(1 - 7/ ^l/4-V™ 2 -™+ 1 /4/2 

r(ai)r(ft) 1 J 

x 2-^1(71 - ai,7i - Pi\1i + 1 - «i - A; ! — «)]• 

Due to the MTBH being asymptotically anti-de Sitter 
the QNM boundary conditions at infinity demands that 
fx — ^ as it — y 1 [37]; [33- From the expression (|2"T)l we 
note that the first term in square brackets vanishes as u — > 
1. The second term vanishes for 1/2 > \J m 2 — rh + 1/4, 
(that is, for 1 > rh). Thus the function fx goes to zero as 
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u — > 1, and therefore if 1 > to the boundary condition at 
infinity docs not impose any restriction on the frequencies, 
that is, there is a continuum of frequencies that satisfy the 
boundary condition at infinity of the QNM. For m > 1, in 
order that fa — > as u — > 1 we must impose the condition 

a\ = — m, or /3i = — m, 7ii = 0,l,2,... (28) 

Therefore for to > 1, from Eqs. (28) we find that the QNF 
of the function fa are equal to 



Ul!=K 



2m 



K — i 2n\ + \ m 2 — rh + 



(29) 



whereas for 1 > to there is a continuum of QNF. 

To calculate the QNF of the function fa we choose the 
quantities C 2 , B 2 , and F 2 as C 2 = iK/2, B 2 = iw/2 + 1/4, 
and F 2 — 1/4 + ^m 2 + to + 1/4/2. A similar method to 
that used for the function fa allows us to find that for all 
to the QNF of the function fa are (n 2 = 0, 1, 2, ... ) 



oj 2 = K — i i^2n 2 + rh 2 + to + 



Cj 2 = -K -i\2n 2 + l + Jm 2 +m+\ . (30) 



From formulas (IT3) we find that the functions R\ and R 2 
are linear combinations of the functions fa and fa, there- 
fore only the QNF that are equal for both functions fa 
and fa will be QNF of the Dirac field in Z?-dimensional 
MTBH. Thus when to < 1, for the function fa we find a 
continuum of QNF, but for the function fa we only find 
the QNF ||3D). Hence for m < 1 the QNF of the Dirac field 
are equal to 



K i / 1 

oj = \ 2n H h to 

L L\ 2 



K i f n 3 „ 



n = 0,l,2, 



(31) 



When to > 1 for the function /i we find the QNF (29) . 
whereas for the function fa we find the QNF ([30) . Af- 
ter some simplifications we find that for to > 1 the QNF 
frequencies of the Dirac field are also determined by the 
expressions (f3"Tj) . Thus in MTBH formulas (|3"Tj) give the 
QNF of the Dirac field for any value of the mass to. In the 
massless limit the QNF ([31]) reduce to 



K 
1 



2n- 



K 



u3 = -T~l 



2n 



(32) 



For QNF (gD and (32) we find that Im (ui) < 0, hence 
these QNM decay in time. Thus the Z?-dimensional MTBH 
is linearly stable against Dirac perturbations. Something 



similar happens for the QNF of the electromagnetic and 
gravitational perturbations [37] , [39] . The stability of the 
MTBH against the gravitational perturbations was shown 
in Refs. gD], [S]. 

As we previously commented, in Refs. |37]-|39) it was 
calculated the QNF of the gravitational, electromagnetic, 
and minimally coupled massless Klein-Gordon perturba- 
tions. The values obtained for the QNF of these fields are 



± 



£ 2i 



(33) 



where the quantity A takes on the values 



D 



1 for the vector type gravitational 
and electromagnetic perturbations, 
\D — 5 1 + 2 for the scalar type gravitational 
and electromagnetic perturbations, 
D + 1 for the tensor type gravitational 

perturbation (D > 5) and minimally 
coupled massless Klein-Gordon field, 

(34) 

the quantity £ depend on the perturbation type and is 
related to the eigenvalues of the Laplacian on the manifold 

s fl _ 2 G2H3S1. 

For the non-minimal coupled to gravity massive Klein- 
Gordon field the QNF are equal to [38] 



(35) 



where 



l eff 



1 — jD (D — 2)/ (4L 2 ) , m denotes the mass 



of the Klein-Gordon field, and 7 is the coupling constant 
between the scalar curvature and the Klein-Gordon field. 
Notice that in Ref. [38] it was chosen a different time pa- 
rameter to that used in the present paper. This fact implies 
that the QNF ([3"5|) have an additional factor of 1/L to the 
QNF reported in Ref. [35]. 

From formulas ([3T1) - ([55"1) we find that for the Dirac field 
the imaginary part of the QNF ([3T]) and ([32]) does not de- 
pend on the spacetime dimension, unlike boson fields the 
imaginary part of their QNF (|33[) and ([35]) shows an ex- 
plicit dependence on the spacetime dimension. Thus for 
the Dirac field the decay time Td = l/|Im(w)| depends on 
the mode number n and it does not depend on the space- 
time dimension. In contrast for the boson fields the decay 
time is inversely proportional to the spacetime dimension, 
thus for a given boson field and fixed mode number, the de- 
cay time decreases as the spacetime dimension increases. 
Hence in D-dimensional MTBH the decay time for the 
Dirac field and the decay time for the boson fields show a 
different behavior when the spacetime dimension changes. 

Furthermore, from formulas ([33]) and ([34) for the mass- 
less boson fields with mode number fixed, and for D > 5 
we find that the tensor type gravitational perturbation 
and minimally coupled massless Klein-Gordon field decay 
faster than vector type and scalar type electromagnetic 
and gravitational perturbations. 
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From QNF (32J) and Q33J) we see that for D > 6 the 
decay time of the massless Dirac field is greater than the 
decay time of the massless boson fields. Thus for D > 6 
the massless boson fields decay faster than massless Dirac 
field. Also for D = 5 the tensor type gravitational per- 
turbation and minimally coupled massless Klein-Gordon 
field decay faster than the other massless boson fields and 
massless Dirac field. For D = 4 we find that the minimally 
coupled massless Klein- Gordon field decays faster than the 
electromagnetic, gravitational, and massless Dirac pertur- 
bations. 

It is convenient to note that for the massive Klein- 
Gordon and Dirac fields the imaginary part of the QNF 
depends on the mass of the field. Taking into account for- 
mulas (|3"Tj) and (|35|) we find that if the mass of the Dirac 
and the minimally coupled Klein-Gordon fields are equal 
and the condition 



mL < 



D 



1 



D 
~2 



(36) 



is satisfied, then the minimally coupled Klein-Gordon field 
decays faster than the Dirac field. 

In MTBH the oscillation frequencies of the boson and 
fermion fields do not depend on the mass of the field. For 
the boson fields the oscillation frequencies are determined 
by the eigenvalues of the Laplace operator on the negative 
curvature manifold £d_2, whereas for the Dirac field the 
oscillation frequencies are determined by the eigenvalues 
of the Dirac operator on £rj_2. 

Thus for a complete determination of the QNF (|3"Tj) for 
the Dirac field moving in MTBH, we need to know the 
eigenvalues of the Dirac operator on the base manifold 
£,d_2 with metric AT? D _ 2 . We expect that the event hori- 
zon of a black hole will be a compact and orientable man- 
ifold [55] . For the MTBH the negative curvature manifold 
£_d-2 usually is a quotient of the form H D ~ 2 /G, where G 
is a freely acting discrete subgroup of the isometry group 
for the (D— 2)-dimensional hyperbolic space H D ~ 2 . There- 
fore for the QNF {ST]) of the MTBH we need to find the 
spectrum of the Dirac operator on a compact spin manifold 
of hyperbolic type. Regarding the spectrum of the Dirac 
operator on hyperbolic manifolds we know the following 
facts. 

In contrast to the Laplace operator, the spectrum of the 
Dirac operator depends on the geometry of the manifold 
and the spin structure, which is a topological object that 
we need to define spinors [B5], [70]. In general, the spin 
structure of a spin manifold is not unique, for example 
the circle S 1 has two spin structures, but note that some 
manifolds do not admit even a spin structure, for example 
the complex projective plane CP 2 [69], [70] . 

The hyperbolic space H has an unique spin structure 
(it is due to that the hyperbolic space is contractible) [7T] . 
It is known that on the hyperbolic space for the Dirac 
operator the discrete spectrum is empty and its continuous 
spectrum is R [69] , [71] . We note that the conventions used 
in Refs. [69], [71] and the present paper are different. In the 



conventions that we use here the eigenvalues of the Dirac 
operator on the hyperbolic space are purely imaginary as 
in Ref. [72] (and therefore K € R), whereas in Refs. [69] . 
[7T] the eigenvalues of the Dirac operator on the hyperbolic 
space are real numbers. 

If the manifold is compact, then general elliptic theory 
asserts that the spectrum of the Dirac operator is discrete 
[69] . [70] . Thus we expect that on the base manifold E.d-2 
of the MTBH the eigenvalues of the Dirac operator are dis- 
crete. Furthermore for a D-dimensional compact manifold 
£ the eigenvalues k of the Dirac operator satisfy the Weyl 
asymptotic law |69] 



lim 



N(k) _ 2[ D / 2 W(£) 



K D (4tt) d / 2 T (f + 1) 



(37) 



where vol(S) is the volume of the D-dimensional manifold 
£ and N(k) is the number of eigenvalues whose modulus 
is < K. 

On a compact symmetric manifold with a homogeneous 
spin structure the square of the Dirac operator J7 2 satisfies 

vr2 



(38) 



where £1 is the Casimir operator of the isometry group and 
1Z is the scalar curvature of the compact symmetric mani- 
fold. Therefore for these manifolds the computation of the 
spectrum for the square of the Dirac operator can be done 
by algebraic methods. Also on these manifolds the spec- 
trum of the Dirac operator is symmetric with respect to 
the origin and the spectrum of the Dirac operator is deter- 
mined by the spectrum of its square. Nevertheless there 
are technical difficulties and the spectrum of the Dirac op- 
erator is explicitly known for a small number of manifolds 
[73]. 

As far as we know for compact hyperbolic manifolds the 
spectrum of the Dirac operator is calculated exactly for 
the manifold £ = PSL 2 (R)/T, where PSL 2 (R) is the pro- 
jective special linear group of R 2 and 1 is a co-compact 
Fuchsian subgroup [73], [74]. The complicated spectrum 
of the Dirac operator on £ = PSL 2 (R)/T appears in Theo- 
rem 2.2.3 of Ref. [73]. Notice that the case relevant to our 
work is when the parameter t of Theorem 2.2.3 is equal 
to 1 and therefore the manifold PSL 2 (M.)/T has negative 
constant sectional curvature!! 

We notice that for the Dirac operator eigenvalue esti- 
mates can be found in several manifolds for which an exact 
calculation of the spectrum is not possible [73] . We believe 
that the following result is relevant for our work. 

In Proposition 2 of Ref. [75] it is asserted that for a 
compact and oriented two dimensional surface £ of genus 



2 For the related case of the so called plane symmetric black hole, 
it is convenient to notice that the spectrum of the Dirac operator on 
the higher dimensional flat tori has been calculated (see Theorem 4.1 
of Ref. [69] and Theorem 2.1.1 of Ref. [73]). We point out that the 
flat tori admits several spin structures and the spectrum of the Dirac 
operator depends on the spin structure [73) . For other examples 
of flat manifolds for which the spectrum of the Dirac operator is 
calculated exactly see Chapter 2 of Ref. |73| . 
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g 7^ 1 there is an eigenvalue K of the Dirac operator that and from the QNF ([321 of the massless Dirac held we hnd 
satishes 



K |< c{g) max{principal curvatures of £}, (39) 



where 



c(g) 



1 if g = 0, 

3 if <? = 2,3, 

4 if 5 > 4. 



(40) 



For g > 2 this result is pertinent for the four-dimensional 
MTBH. We do not know similar estimates for the eigenval- 
ues of the Dirac operator on higher dimensional compact 
hyperbolic manifolds. 

From these comments it is deduced that in the mathe- 
matical literature we do not find many calculations on the 
eigenvalues of the Dirac operator on compact hyperbolic 
manifolds and we believe that the computation of these 
quantities is a challenging mathematical problem. 

In Ref. [5] it was shown that the momentum space poles 
of the retarded correlation functions in the dual conformal 
field theory and the QNF of the three-dimensional BTZ 
black hole are identical. Calculating whether something 
similar happens for the QNF of the D-dimensional MTBH 
is an interesting problem. 



ttT h 



= 1 



and 



Hull 
ttT h 



(44) 



We see that the first expression in formulas (|44j) saturates 
the inequality (j42|) and that the second expression does 
not satisfy the previously mentioned inequality. 

Furthermore, from QNF ([33]) of the massless boson fields 
we obtain that 



ttT h 



(45) 



Hence, taking into account the values of the quantity A 
given in formula ([51]). for D > 4 we find that in MTBH 
the fundamental QNF of the massless bosons do not satisfy 
the inequality (|4"2"]> . 

Thus we find that in D-dimensional MTBH the funda- 
mental QNF of the massless boson and Dirac fields do not 
satisfy inequality (|42|) . We expect that inequality (|42|) be 
satisfied in MTBH [58] owing to Hawking's temperature 
of the MTBH is of the same order of magnitude as the re- 
ciprocal of the characteristic length (L) of the spacetime. 
According to Hod the TTT bound fill) is universal and we 
do not know the cause of its failure for the fundamental 
QNF of the MTBH. 



3 Hod's bound 



4 Effective potentials 



Taking into account quantum information theory and ther- 
modynamic concepts, in Ref. [58] Hod found a bound on 
the relaxation time r of a perturbed thermodynamic sys- 
tem. This bound is 



h 

7rT' 



(41) 



where T m i n stands for the minimum relaxation time and 
T denotes the temperature of the thermodynamic system. 
This bound is called "TTT bound" (time times tempera- 
ture bound) by Hod in Ref. [58 . 

In Ref. [58] it was shown that strong self-gravity sys- 
tems, as the black holes, are the appropriate systems to 
test the TTT bound (01]). For a black hole the TTT bound 
states that at least for the fundamental QNF the following 
inequality is satisfied [55] 



< 1, 



H 



(42) 



where wj is the absolute value of the imaginary part of the 
fundamental QNF and Th is Hawking's temperature of the 
black hole (see Refs. [58], [76]-[79] for more details). The 
fundamental QNM is the least damped mode of the black 
hole and it determines its relaxation time scale PQ-[3J. 
The Hawking temperature of the MTBH is equal to [50]- 

T « - £z> (43) 



Following the method of Chandrasckhar's book [55J, we 
take for the Dirac field a harmonic time dependence as in 
formula ([9]) to transform Eqs. (JSJ into the pair of decoupled 
Schrodinger type equations 



d 2 Z ± 
df* 2 



■ to Z± = V± Z± , 



(46) 



where 



Z ± = C l0 / 2 R 1 ±e-^l 2 R 2l 
mz 

= arctan^^, 
K 

, dW 

v± = w 2 ±—, 



(47) 



3/2 



w 



Vz^l(k 2 + (mz) 2 ) 

z(ir 2 + (m^) 2 ) + ggz(z 2 -i) : 



K = K/L, we define r* by 



df* 
dr* 



mK 



2loL K 2 + (r 



(48) 



and, as in Sec. [2] r* denotes the tortoise coordinate of the 
.D-dimensional MTBH (see formula (f2"T]l ), 

From formulas (|46|) and (|47]) we see that the effective 
potentials V± are complicated functions of the different 
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parameters. Nevertheless in the massless limit we find that 
the formulas for W and V± reduce to 



W = K- = -A'scch(r*/L), 

V ± = f ± (k/L)^/L)_ 
cosh {r*/L) cosh (r*/L) 

Thus for the massless Dirac field the effective potentials 
igS} are of Morse type [50]. In Ref. [35] it was shown 
that in £>-dimensional MTBH the effective potentials of 
the Schrddinger differential equations for the massless bo- 
son fields are of Poschl- Teller type. We note that for many 
of the spacetimes for which we exactly calculate their QNF 
the effective potentials of Schrodinger type equations are 
of Poschl- Teller o Morse type (see for example Table 1 in 
Ref. [23]). 

5 Discussion 

For the D-dimensional MTBH in Sec. [2] we found that the 
real part of the QNF is determined by the eigenvalues of 
the Laplace operator (boson fields) or the Dirac opera- 
tor (Dirac field) on the negative curvature manifold Sd-2- 
Nevertheless, up to our knowledge there are no many cal- 
culations of the spectrum of the Dirac operator on com- 
pact spin manifolds of hyperbolic type. We believe that 
this mathematical problem deserves detailed study. Fur- 
thermore we notice that the imaginary part of QNF ([31]) 
is independent of the eigenvalues of the Dirac operator on 
T,d-2- This fact allows us to discuss some phenomena (see 
Sees. [5] and [3J even if we do not know explicitly the value 
of the eigenvalues of the Dirac operator on the manifold 

For the massless boson and Dirac fields the imaginary 
part of the QNF shows a different dependence on the space- 
time dimension. For the boson fields the decay time de- 
pends on the spacetime dimension whereas for the Dirac 
field it is independent of the spacetime dimension. Also we 
point out that for D > 6 the massless boson fields decay 
faster than the massless Dirac field. 

In MTBH the QNF of the Klein-Gordon, gravitational, 
electromagnetic, and Dirac perturbations have been calcu- 
lated (see Refs. [37]-[39] and Sec.[2]of this paper). Never- 
theless as far as we know the QNF of the Rarita-Schwinger 
field have not been computed. We believe that the calcu- 
lation of the QNF for this field is an interesting problem. 

According to Hod the TTT bound of formula (gT]) is 
universal [55], [ZS] _ [ZS], but we found in Sec. [3] that for 
the fundamental QNF of the MTBH the inequality flU) is 
not satisfied (see formulas ([44]) and (J45J) ) . We believe that 
this puzzling result deserves detailed study. 

For the Z3-dimensional MTBH, from our results and 
those already published we obtain that the real part of 
the QNF depends on the eigenvalues of the Laplace or 
Dirac operators on the negative curvature manifold T,o-2- 
These values can be different for distinct fields, also for a 



fixed field these eigenvalues may depend on the mode of 
the field. Thus the asymptotic limit of the real part of the 
QNF for the D-dimensional MTBH depend on the phys- 
ical parameters of the black hole and the field (and the 
mode of the field) . 

An interesting proposal is the so called Hod's conjecture 
[6], it states that in the semiclassical limit the area quan- 
tum of a event horizon can be calibrated with the asymp- 
totic value of the real part of the QNF. The facts men- 
tioned in the previous paragraph imply that Hod's con- 
jecture is not valid for the D-dimensional MTBH (as for 
the D-dimensional de Sitter spacetime [81]), since in this 
conjecture we must assume that the real part of the QNF 
depends only on the physical parameters of the black hole 
Eh 0; but it does not happen in £>-dimensional MTBH. 
Thus we think that for the D-dimensional MTBH we must 
investigate whether the recent proposal of Maggiore [7] can 
be used to determine the area quantum of its event hori- 
zon. Work along this line is in progress. 

Finally we notice that formulas (J3TJ) also give the QNF 
of the Dirac field propagating in three-dimensional static 
BTZ black hole with mass M = 1. The QNF of the Dirac 
field evolving in the static BTZ were previously calculated 
in Refs. [5], p]. 
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